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ABSTRACT:  Inelastic  scattering  of  atoms  of  moderate  energies  (say  <  5  eV)  by  solid 
surfaces  is  almost  entirely  due  to  energy  exchange  with  lattice  vibrations.  It  can  give  valu¬ 
able  information  about  the  atom-surface  interaction  potential  and  the  vibrational  dynamics 
at  surfaces.  Theoretically  this  process  represents  a  challenging  many-body  problem,  call¬ 
ing  for  suitable  approximation  methods.  In  this  paper  I  shall  outline  work  currently  in 
progress  in  our  group  (K.  Burke,  L.  D.  Chang,  W.  Kohn).  (l.)  We  have  solved  a  simple 
model  problem  in  which  the  normal  modes  of  the  lattice  are  schematized  by  a  single  one¬ 
dimensional  harmonic  oscillator,  initially  in  its  groundstate  (T  =  0).  The  classical  solution 
gives  a  unique  energy  loss.  We  have  calculated  the  leading  quantum  correction  and  find 

a  Gaussian  final  energy  distribution  whose  width  is  proportional  to  h.2.  Our  exact  results 
are  in  general  different  from  the  so-called  trajectory  approximation.  (2.)  We  are  about 
to  propose  a  new  type  of  atom-surface  scattering  experiment,  which  will  provide  a  direct 
measure  to  the  quantum  corrections  to  classical  scattering. 
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ABSTRACT:  Inelastic  scattering  of  atoms  of  moderate  energies  (say  <5  eV) 
by  solid  surfaces  is  almost  entirely  due  to  energy  exchange  with  lattice  vibrations. 
It  can  give  valuable  information  about  the  atom-surface  interaction  potential  and 
the  vibrational  dynamics  at  surfaces.  Theoretically  this  process  represents  a  chal¬ 
lenging  many-body  problem,  calling  for  suitable  approximation  methods.  In  this 
paper  I  shall  outline  work  currently  in  progress  in  our  group  (K.  Burke,  L.  D. 
Chang,  W.  Kohn).  (1.)  We  have  solved  a  simple  model  problem  in  which  the 
normal  modes  of  the  lattice  are  schematized  by  a  single  one-dimensional  harmonic 
oscillator,  initially  in  its  groundstate  (T  =  0).  The  classical  solution  gives  a  unique 
energy  loss.  We  have  calculated  the  leading  quantum  correction  and  find  a  Gaus¬ 
sian  final  energy  distribution  whose  width  is  proportional  to  /12 .  Our  exact  results 
are  in  general  different  from  the  so-called  trajectory  approximation.  (2.)  We  are 
about  to  propose  a  new  type  of  atom-surface  scattering  experiment,  which  will 
provide  a  direct  measure  of  the  quantum  corrections  to  classical  scattering. 
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Preface: 


My  friendship  with  Howard  Reiss  goes  back  over  a  little  more  than  one-half  of  our  lives, 
to  the  summer  of  1953.  I  had  come  to  Bell  Labs  on  the  first  of  about  a  dozen  annual  summer 
visits  and  found  Howard  in  masterly  command  of  the  theory  of  solid  solutions  of  Group  in 
and  Group  V  atoms  (acceptors  and  donors)  and  their  complexes  in  Si  and  Ge.  This  was,  of 
course,  during  the  heady  years  of  the  spectacular  transistor  revolution,  and  Howards  work 
was  at  the  theoretical  center  of  the  preparation  of  the  essential  n-  and  p-doped  samples 
of  Si  and  Ge.  As  a  neophyte  in  solid  state  physics  I  was  enormously  impressed  (not  to 
say  intimidated)  by  the  breadth  of  understanding  and  variety  of  techniques  which  Howard 
had  at  his  command. 

In  the  intervening  years  we  both  moved-along  quite  different  paths-into  faculty  po¬ 
sitions  at  the  University  of  California,  and  there  have  been  many  occasions  when  I  have 
enjoyed,  and  frequently  profited  from,  Howard’s  extraordinary  qualities  as  a  scientist.  His 
great  originality  and  outstanding  ability  to  solve,  somehow,  almost  any  reasonable  problem 
in  physical  or  solid-state  chemistry  are  truly  astonishing. 

Howard  is  equally  at  home  in  classical  and  quantum  theory.  I  hope,  therefore,  that  the 
following  contribution,  which  is  on  the  borderline  between  classical  and  quantum  physics 
(or  chemistry)  will  be  to  his  liking. 

I.  Introduction 

In  many  branches  of  quantum  physics  a  semi-classical  description  has  been  the  most 
appropriate  and  most  instructive.  Two  examples  are:  bound  states  in  atoms  and  molecules, 
of  high  quantum  numbers  n;  and  the  description  of  crystal  electrons  in  a  weak  magnetic 
field  H,  where  the  large  parameter  is  the  label  n  of  the  Landau  level  (or,  in  complicated 
geometries,  something  similar).  Semiclassical  treatments  of  scattering  of  a  particle  by 
a  fixed  external  potential  are  also  long  familiar  and  governed  by  the  BWK  equation  or 
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semiciassical  path  integral  methods. 

The  problem  of  inelastic  scattering  of  an  atom  by  a  surface  is  conceptually  and  calcula- 
tionally  much  more  difficult.  If  the  velocity  of  the  incident  atom  (“incidon”)  is  sufficiently 
smaller  than  typical  electronic  velocities-and  this  is  true  up  to  ~lfc  tV  for  H  and  ~100  k  tV 
for  heavy  atoms-the  Born  Oppenheimer  approximation  is  valid  and  the  system  can  be  de¬ 
scribed  by  a  Hamiltonian  of  the  form 

„2  '  "2 


(i) 


where  (m,p,  r)  describe  the  incidon,  and  (M,  p},  ry)  the  N  lattice  atoms.  Thus  electronic 
degrees  of  freedom  are  completely  represented  by  the  interaction  potentials  W  and  V.  We 
shall  limit  ourselves  to  this  Born-Oppenheimer  regime. 

Provided  that  the  surface  temperature,  T,  is  much  greater  than  the  Debye  temperature 
&D  and  the  incident  atom  has  a  DeBroglie  wavelength  much  shorter  than  a  characteristic 
length  of  the  surface  (say  <<10-8  cm)  a  fully  classical  description  evidently  is  in  order. 
If  these  conditions  are  almost,  but  not  quite,  satisfied  a  semiciassical  approach  is  clearly 
appropriate. 

But  what  if,  initially,  the  surface  temperature  Tis  smaller  than  the  Debye  temperature 
so  that  at  least  the  high-lying  modes  are  only  weakly  excited  (n  =  0  or  1)  before  the 
collision;  or,  a  more  extreme  case,  T  =  0,  so  that  initially  al]  modes  are  in  their  n  —  0 
ground-states?  Is  the  classical  (normally  n>>l)  limit  at  all  am  appropriate  starting  point? 
And  what  are  the  appropriate  conditions,  if  any,  on  the  masses  m  and  M,  the  incident 
velocity  v,  and  other  parameters  such  as  interaction  strengths  and  ranges? 

A  good  deal  of  progress  has  been  made  in  recent  yearns  on  semiciassical  inelastic  scatter¬ 
ing  theory.  It  is  of  great  interest  not  only  for  atom-surface  scattering  but  also  for  scattering 
of  molecules  by  surfaces,  accompanied  by  internal  molecular  excitations,  and,  of  course, 
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related  to  chemical  reactions  between  gaseous  molecules.  To  my  mind,  the  above  and  other 
questions  have  not  yet  been  fully  answered  and  considerable  work  remains  to  be  done. 

1.  A  Simple  Model  and  a  Check  on  the  Trajectory  Approximation. 

In  recent  work  by  our  group  we  have  critically  examined  a  simple  model  in  which  the 
atoms  of  the  surface  are  schematized  as  a  single  one-dimensional  harmonic  oscillator.  Our 


Hamiltonian  is 


H  =  L  +  {m  +  lxi)^~{x~x)la 


By  a  simple  change  of  variable,  x—*x  +  b ,  we  have  set,  without  loss  of  generality,  the 
interaction  amplitude  equal  to  unity.  Similar  models  have  been  studied  extensively  in  the 
past  from  various  points  of  view.1  Our  objectives  have  been  two-fold.  (1.)  We  start  with 
the  harmonic  oscillator  in  the  ground  state  and  calculate  the  energy  loss  spectrum  of  the 
incidon  exactly  up  to  and  including  the  leading  quantum  mechanical  correction  (0(h1/2)). 
(2.)  We  compare  these  results  with  those  of  the  so-called  trajectory  approximation2  (which 
has  semiclassical  features)  to  determine  under  which  conditions,  if  any,  that  approximation 
is  valid. 

In  our  model  there  are  2  dimensionless  classical  parameters:  and  p 2=wr, 

where  ui  is  the  frequency  of  the  oscillator  (w  =  /C1/2Af-1/2)  and  r  is  the  time  taken  by 
the  incidon,  of  energy  E,  to  traverse  the  distance  a  (r  =  a  m1//2(2J£)-1/2). 

Our  method  of  solution  is  to  develop  a  singular  limiting  solution  of  the  Heller  wavepacket 
equations.3  Heller’s  method  starts  with  the  oscillator  in  its  Gaussian  groundstate  and  the 
incidon  described  by  a  distant,  incoming  Gaussian  wave  packet  of  the  form 

0(x,t)  =  A I  dpe-(P-Po)3/^e-i(p(*-*o)+W)  (3) 

Here  pQ  and  x0  are  the  initial  mean  momentum  and  position;  £’(p)=p2/2m;(ha)1/2  is  the 
momentum  width  of  the  incident  wave-packet;  and  A  is  the  appropriate  normalization 
constant. 
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The  following  heuristic  argument  suggests  (correctly,  as  we  have  confirmed)  that  the 
momentum  (and  energy)  distribution  of  the  reflected  particle  will  have  a  width  och1/2  for 
h—+ 0:  The  target  oscillator  has  a  zero-point  momentum  distribution  of  width  AP~(hu>M) V2 
which,  using  classical  mechanics,  leads  to  a  distribution  of  reflection  momenta  of  Ap-APa/i1/2. 
To  distinguish  this  real  quantum  spread  from  the  artificial  momentum  spread  of  the  inci¬ 
dent  wave-packet  we  must  make  the  latter  relatively  negligible,  i.e., 

[ha)^^«f{pm 2)(huM)1/2  or 

a«(Mw)1/2  /(Ml, M2)  (4) 

where  /  is  some,  as  yet  unknown,  function  of  the  dimensional  parameters  m  and  p. 2- 
Heller’s  equations  become  singular  in  the  limit  a— >0  and  had  to  be  solved  with  great  care 
by  an  expansion  in  powers  of  a.  The  result  is  that,  up  to  order  h1/2,  the  center  of  the 
reflected  energy  distribution  is  equal  to  the  classical  reflection  energy;  and  that  the  width 
is  given  by 

&E  =  F(mi,m2)(M1/2>  (5) 

where  F(p\,p 2)  is  a  function  of  the  dimensionless  variables  p\  and  P2>  which  we  have 
calculated  for  various  values  of  p\  and  P2  from  our  expansion  of  Heller’s  equations. 

We  have  compared  these  results  with  those  of  the  so-called  trajectory  approximation 
and,  for  any  finite  values  of  p\  and  P2»  find  quantitative  disagreements,  from  a  few  percent 
to  factors  of  the  order  of  2.  Thus  we  conclude  that,  while  the  trajectory  approximation 
provides  a  useful  orientation  concerning  the  magnitude  of  quantum  effects,  it  is  not  quan¬ 
titatively  correct. 

2.  Specific  Quantum  Effects  in  Atom  Surface  Scattering  Experiments. 

If  an  atom  is  incident  on  a  perfect  surface  at  temperature  T  with  a  deflnite  momen¬ 
tum  pi,  its  reflected  momentum  pj  will  have  a  spread  in  both  direction  and  magnitude. 


Classically  this  spread  is  due  to  the  different  points  of  impact  and  to  the  different  ini¬ 
tial  states  of  the  target  in  the  temperature  ensemble.  In  addition,  quantum  mechanically, 
even  for  a  given  initial  state  of  the  target  there  are  many  possible  final  target  states  and 
corresponding  final  momenta  pj.  To  isolate  this  last,  specific  quantum-mechanical  effect 
we  are  proposing  the  following  experiment:  The  target  should  be  at  temperature  0,  so 
that  variations  of  the  initial  target  state  are  negligible.  Secondly  the  reflected  incidon  is  to 
be  observed  in  a  definite  direction,  Qj.  Classically  one  expects  exactly  2  sharp  reflection 
energies,  E\(0  f),E2{0f)  corresponding  to  precisely  2  inequivalent  impact  points  leading 
to  scattering  in  the  chosen  direction  0 j.  (This  effect  is  due  to  the  warping  of  the  surface.) 
If  the  quantum  corrections  are  small  each  of  these  2  sharp  peaks  will  become  a  Gaussian 
line  of  finite  width  ah1/2. 

We  are  currently  studying  which  surfaces,  incident  atom  energies  and  directions  of 
incidence  promise  to  show  this  predicted  quantum  phenomenon  most  clearly. 

Concluding  Remarks. 

Full  accounts  of  this  work  will  be  published  in  the  near  future.  I  would  like  to  acknowl¬ 
edge  the  collaboration  of  L.  D.  Chang  in  the  work  reported  in  Section  1  and  discussions 
with  K.  Burke  concerning  Section  2.  H.  Metiu  has  introduced  us  to  the  Gaussian  wave- 
packet  techniques  and  on  many  occasions  given  us  the  benefit  of  his  experience.  D.  N. 
Newns  has  acquainted  us  with  the  trajectory  method.  Finally  I  acknowledge  with  thanks 
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